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ULTRASOUND MODULATED BIOLUMINESCENCE 
TOMOGRAPHY AND CONTROLLABILITY OF THE RADIATIVE 

TRANSPORT EQUATION 


GUILLAUME BAL, FRANCIS J. CHUNG, AND JOHN C. SCHOTLAND 

Abstract. We propose a method to reconstruct the density of an optical source 
in a highly scattering medium from ultrasound-modulated optical measurements. 
Our approach is based on the solution to a hybrid inverse source problem for the 
radiative transport equation (RTE). A controllability result for the RTE plays an 
essential role in the analysis. 


1. Introduction 

This paper is concerned with the problem of reconstructing an optical source in a 
highly-scattering medium. The primary application is to bioluminescence imaging, 
in which the cells of a model organism are tagged with a light-emitting molecular 
probe [m EH]. The goal is to recover the spatial distribution of the labelled cells 
from measurements of the emitted light. This is a classical inverse source problem 
for the radiative transport equation (RTE), which may be formulated as follows. 

Let A be a bounded domain in M” with smooth boundary dX, for dimension n >2. 
The specihc intensity u obeys the RTE 

9-Vu + (t{x)u— f k{x,9,9')u{x,9') d9' = S{x) . (1.1) 

J S’^-^ 

Here u{x, 9) is the intensity of light at the point x ^ X traveling in the direction 
9 G S^~^, S is the source (which is taken to be isotropic) and a is the attenuation 
coefficient, which is taken to be nonnegative. The scattering kernel k is nonnegative 
and obeys the reciprocity relation k{x, 9, 9') = k{x, —9', —9) and the normalization 
condition J k{x,9,9')d9' = 1 for all 9 G We will also assume that k is invari¬ 

ant under rotations. That is, k{x,9,9') = k{x,9 ■ 9'), which holds for statistically 
homogeneous random media. We dehne the subsets r± of dX x S"'~^ by 

r± = {(x, 9)edX X : ±9 ■ n{x) > 0} , (1.2) 

where n{x) is the outward unit normal vector at x. The specific intensity obeys 
the boundary condition u{x, 9) = 0 for [x, 9) G r_. Thus no light enters the domain 
except from the source. The inverse source problem is to reconstruct S from boundary 
measurements of the outgoing specihc intensity M|r+; assuming that a and k are 
known. This problem has been considered in [12] , where it is shown that it is possible 
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to uniquely recover S, provided that k is sufficiently small in a suitable norm. See 
also [231 EH] for related work. In [3D] , it was shown that the smallness condition can be 
replaced by the assumption that a and k belong to C^{X). A corresponding stability 
estimate was also derived. However, the case of continuous coefficients remains open. 

Ultrasound modulated bioluminescence tomography (UMBLT) is a recently pro¬ 
posed imaging modality in which an acoustic wave is used to spatially modulate the 
optical source, while boundary measurements of the optical held are performed [m 
[20] . In this manner, the above inverse source problem is converted into a so-called 
hybrid inverse problem, where an external held is used to control the material prop¬ 
erties of a medium of interest, which is then probed by a second held. See [SHHISIISI 
OIHIISI HDl EDI EDI Ell EDI EH ESI ESI EZ] for examples of hybrid inverse problems in 
other physical settings. 

The inverse problem of UMBLT was hrst studied in the special case of the dihusion 
approximation (DA) to the RTE [11]. It was shown that it is possible to uniquely 
reconstruct the source with Lipschitz stability. However, the DA breaks down in 
strongly-absorbing or weakly-scattering media, near boundaries and on small length 
scales. In this paper, we consider the inverse problem of UMBLT in the transport 
regime. Once again, we hnd that the source can be recovered with Lipschitz stability 
by a constructive procedure. 

The forward problem of UMBLT is formulated as follows. Suppose that a standing 
acoustic pressure wave with modulation amplitude of the form cos{q-x + (p) is incident 
on a highly-scattering medium, where q is the wave vector and ip is the phase of the 
wave. Following mm, we hnd that the the coefficients a and k and the source S 
are modulated thus forming new coefficients Ug, ke and S^, which are given by 

(Js{x) = (1-f £ cos(g • X-|- 9 ?))(T(a;) , (1.3) 

k^{x) = {1 + e cos{q ■ X + ip))k{x) , (1.4) 

Se{x) = {1 + e cos{q ■ X + ip))S{x) , (1.5) 

where 0 < e 1 is the dimensionless amplitude of the acoustic wave. The RTE thus 
becomes 

9 ■ Vue + ae{x)u — f k^{x,9,9')ue{x,9') d9' = Se{x) , (1.6) 

where the dependence of the specihc intensity u on e has been made explicit. The 
inverse problem now consists of deducing S from boundary measurements of Ue when 
q and (p are varied. 

Throughout this paper we will make the following assumptions on the domain X 
and the coefficients a and k that appear in the RTE. First, we will assume that X is 
a bounded subset of M” with diameter r. Next, we will assume that a, fc, and S are 
continuous. Finally, in order to ensure solvability of the RTE, let 

(1.7) 

and assume that one of the following inequalities holds: 


P = 




k{x,e,e')de' 


a — p > a, 


(1.8) 
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for some positive constant a, or 

rp < 1. (1.9) 

Under the above conditions, the equation fll.hp with boundary condition Me|r_ = 0 
has a unique continuous solution (see Section 2 for details), and we can dehne 
V^) : M"- X {0, |} —)■ C(r+) to be the map given by 

-“o)|r+ , (1.10) 

if £ 7 ^ 0. (Recall that is defined in part by the choice of q and p.) We will also need 
to make separate use of the measurements when e = 0, that is without ultrasound 
modulation. Therefore we will define 

As = ^o|r+- (1-11) 

We have the following theorem. 

Theorem 1.1. Suppose X, a and S are as above and that k is invariant under 

rotations. Then the map S eA {A^ : e > 0} zs injective. Moreover, we have the 

stability estimate 

UlIRi - S'2 ||l°°(a) < ll^si “ ^5211^^+) + 11^5^ - A5j|Li(Rnx{o,f},c(r+)) + 0(e) (1.12) 

where C depends only on X, k, and a. 

The proof of Theorem 11.11 is constructive. That is, we will also give an algorithm 
by which S can be reconstructed from knowledge of A|. 

Remark 1.2. Note that as compared to the result o/pSU], the above stability estimate 
does not depend on derivatives of the coefficients a and k. 

A key ingredient in the proof of Theorem 11.11 is the following controllability result 
for the RTE, which is of interest in its own right. 

Theorem 1.3. Suppose X, a, and k are as given above. Then for any point xq G X 
and any continuous function h onU’{S'^~^), there is a function g G Lp{S"'~^, L°°{dX)) 
such that the boundary value problem 

9-Xv + av = f k{x,9,9')v{x,9') d9' 

Js^-^ 

w|r_ = 5'|r_ 

has a unique solution v G {S^~^, L°°{X)) which is continuous in a neighbourhood 
of Xq, and satisfies the property that v{xo,9) = h{9), for all 9 G S"‘~^. Moreover, for 
any 1 <p < oo, 

\\v\\lp(S"-^,l°°(x)) + ||5'||LP(S"-i,L°°(av)) < C\\h\\LP(s"-i). 

where C depends only on X, a, and k. 


(1.13) 
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The proof of Theorem 11.31 is also constructive. Note that other control theory 
results for the time-dependent RTE have had applications to inverse problems as 
well. See for example BE]. 

The remainder of this paper is organized as follows. In Section 2, we prove a 
regularity result for the RTE which will be needed in the proofs of Theorem 11.11 and 
Theorem 11.31 We then prove Theorem 11.11 in Section 3, using Theorem 11.31 Finally, 
we will prove Theorem 11.31 in Section 4. 


2. Regularity 

In order to prove Theorem 11.11 and Theorem 11.31 we need the following regularity 
result for the RTE. 


Theorem 2.1. Let /_ G (^(r.). Under the conditions on X, a, k, and S given in 
Section 1, the equation fll.ip has a unique continuous solution u with the boundary 
condition 

w|r_ = /- (2.1) 

Moreover, if fll.81) holds, then 

||w||lp(S'‘-i,C(A)) < —((p + a)||/-||LP(S'‘-i,C(aA)) + ||*S'||LP(5n-yc'(X)))- (2.2) 

If instead fll.91) holds, then 

||w|1lp(5’"-i,C(A)) < ^ _ ^^ (Il/-||LP(S’”-I,c(ax)) + l||5'||lp(S"-i,C(X)))- (2.3) 

In both cases, the estimate given is valid for 1 < p < oo. 


The proof of Theorem 12.11 is close to that of similar results in [12]; see also 


Let 

T±{x, 9) = min{t >0\x±t9E dX] 

(2.4) 

and 

Aiu = au , 

(2,5) 


A 2 U = — f k{x,9,9')u{x,9')d9' , 

Js^-^ 

(2.6) 


Au = Aiu + A 2 U , 

(2,7) 


Tqu = 9 ■ Xu , 

(2.8) 


Tiu = To -|- y4i. 

(2,9) 

In addition, let 




B{t, X, 9) = exp f — J a{x — s9, 9)ds\ . 

(2.10) 


Note that B is continuous in each variable and \B{t,x,6)\ < 1. Moreover, if fll.81) 
holds, then we have the improved estimate 

\B{t,x,9)\<e-^^P+'^l ( 2 . 11 ) 

Finally, if /_ G C(r_), we dehne J/_ to be the function on X x S"‘~^ dehned by 

Jf- = B{T-{x,9),x,9)f-{x - t-{x,9)9,9). (2.12) 
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Note that J/_ is continuous, TiJ/_ = 0, and J/_|r„ = /-• Moreover, for any fixed 

6 e 

||-^/-||c(A) < ||/-||c(9X), (2-13) 

and so 


\\jf-\\LP(S”-^,C{X)) < ||/-||lp(S'"-i,C(9A))- 


(2.14) 


Proof of Theorem \2.1i Suppose u satisfies fll.ll) with the boundary condition M|r_ = 
/_. Then on X x we have in the above notation that 


(Ti + A2)u = S. 

Thus for any t G [0, T-{x, 9)], we have 

B{t, X, 0)[(Ti + A 2 )u]{x — to, 9) = B{t, X, 9)S{x — tO). 

Now let 

Tf^f{x,9)= / B(t,x,9)f{x — t9,9)dt. 

Jo 

Then T^^f is continuous if / is. Moreover, for any hxed 9, 


f\\c{x) < p||/||c(A), 


and if we assume fll.Sp . then 


rT-(x,e) 


Ti f\\c{x) < ||/||c(A) / \B{t,x,9)\dt 
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< 


C{X) 




< (P + a) ll/llc(A)- 


Therefore it follows that 


(2.15) 

(2.16) 

(2.17) 

(2.18) 

(2.19) 

( 2 . 20 ) 
( 2 . 21 ) 


11^1 V||lp(S'*-i,C(X)) < p||/||lp(S'*-i,C(X))5 (2.22) 

and if fll.Sp holds, then 

11^1 V||lp(S"-i,C(A)) < (p + «) ^||/||lp(S"-i,C(A))- (2.23) 

If we now integrate fl2.16p in t from 0 to T-{x,9), then we obtain after some 
calculation the result 

(/ + Tf^A2)u = Jf. - T^^S. (2.24) 

Therefore, if u solves (II.ip with boundary condition M|r_ = /-, then u solves (I2.24p . 
Moreover, if we apply Ti to both sides of (12.241) . then we get (II.ip . and since T-{x, 9) = 
0 on r_, it follows that on r_, (12.241) reads u = f-. Therefore u solves (11.11) with the 
boundary condition M|r_ = /_ if and only if u solves (12.241) . 

Now if K = Tf^A 2 , then (I2.24p reads 

(J + K)u = Jf. - Tf^S. (2.25) 


Moreover, if (11.81) holds, then 
\\Ku\\lp(^S^-i^C{X)) < 


P 


p + a 


\u 


||lp(S"-i,C(A)) < ||m||lp(S"-i,C(A)), 


(2.26) 
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and if (ll.9j) holds, then 

WKuWlp^^s-^-i fi(x)) < 'TP\\u\\lp{s^-^,c{x)) < \\u\\lp{S'^-^,c{x))- (2.27) 

In either case, (/ + K)u has a unique inverse which can be written using Neumann 
series. The series 

u = {I-K + K‘^-K^ + -- ■){Jf_ - T^^S) (2.28) 

converges uniformly in {x, 6), since we can set p = oo, and hence yields a continuous 
function u, which must satisfy the inequality fl2.2p in the case that fll.Sp holds, and 
fl2.3p in the case that fll.Qp holds. Thus u solves fl2.24p and thus solves fll.ip . 

□ 

As a remark, note that if boundary source /_ G L°°{dX)) but is not con¬ 
tinuous, then the same argument as above will give a solution u G L°°{X)). 


3. Inverse Problem 


We recall that the inverse problem consists of deducing the source S from boundary 
measurements of the specihc intensity u^, assuming the coefficients a and k are known. 
Here we also assume that k is invariant under rotations. To proceed, we suppose that 
n is a continuous solution to the adjoint equation 


— 9-Xv + av= f k{x,9 ■ 9')v{x,9') d9' , (3.1) 

Js”-^ 

obeying the boundary condition n|r+ prescribed. Since k and a are known, we will 
suppose that v is known. We now make use of the integration by parts identity 


v9 ■ Xur dx 


Ue9 ■ Vv dx + 


UeV n ■ 9 dx 


(3.2) 


Jx Jx Jax 

where n is the outward normal vector to dX. If v satishes fl3.ll) and Ue satishes fll.bp . 
then fl3.2l) becomes 


— I VaeUe dx + I v{x, 9) / ke{x, 9 ■ 9')Ue{x, 9') d9' dx + 


> X 


'X 


>x 


Ui.avdx+ / Ui;{x,9) / k(x,9 ■ 9')v{x,9') d9'dx 

Jx 


UpV n ■ 9 dx. 


fax 


Rearranging and integrating over 9 gives 

/ / UpV n ■ 9 dxd9 


lax 


(cr — (Jp)upV dx d9 + 


vSp dx d9 


lx 


lx 


+ 



’X 


X Js^-^ Js^-^ 



X 


v{x, 9)kp{x, 9 ■ 9')us{x, 9') d9' d9 dx 
Up{x, 9)k{x, 9 ■ 9')v{x, 9') d9' d9 dx. 


vSp dx{3.3) 

(3.4) 

(3.5) 

(3.6) 

(3.7) 

(3.8) 

(3.9) 
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We can now write the above as 



u^v n ■ 6 dxdO 


lax 


[cr - a. 


)ueV dx dO + 


vS^ dx dO 


Jx 


lx 


+ 




k){x, 9 ■ 9')v{x, 9)u^{x, 6') dO' dO dx. 

(3.10) 


Since we can measure u^\gx-, and v is assumed to be known, the left hand side of the 
above equation is a known quantity. Now if £ = 0, then writing u = Uq, we have 



uvn ■ 9 dxdO 


lax 



vS dx d9. 


(3.11) 


This fails to take advantage of the ultrasound modulation, though. If we dehne ul 
by 

eul = Us — u, (3-12) 

then ul satishes the equation 


s9 ■ Vul + sa^u] = e 


ks{x,9,9')ul{x,9')d9' 


(3.13) 


+Ss — S+ / {ks — k){x,9 ■ 9')u{x,9') d9'+ {a — as)u. (3.14) 

J s-^-i 


Since the whole second line above is 0(e), it follows from Theorem 12.11 that ||Mg|| is 
0(1). Therefore 


ulvn ■ 9 dxd9 


lax 


JX 


cos(g • X + (p)auv dx d9 + 


(3.15) 

V cos{q ■ X + (p)S dx d9 (3.16) 


Jx 

cos{q ■ X + ip)k{x, 9 ■ 9')v{x, 9)u{x, 9') d9' d9 dx + 0(e). (3.17) 


Jx JS'^-i JS'^-i 

Here the above left hand side is still known. If we now vary q and ip and take e —?■ 0, 
we can obtain the Fourier transform of the quantity H^, dehned by 


H,{x) 


auvd9+ 


vSd9 


(3.18) 


+ 




k{x, 9 ■ 9')v{x, 9)u{x, 9') d9' d9. 


(3.19) 


We note that the above expression for has two unknowns, u and S. However, 
equation fll.ip also relates u and S. If we use fll.ip to substitute for S in the expression 
for Hy, nearly everything cancels and we hnd that 


Hy{x) 



v{x, 9)9 ■ Vui^x, 9) d9. 


(3.20) 
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Remark 3.1. We note that H^^x) is a so-called internal functional; it is known for 
every point x & X. Such internal functionals, which are determined from boundary 
measurements, play the role of internal measurements and are a generic feature of 
hybrid inverse problems. 

The inverse problem now consists of recovering S from Suppose Xq E X. 
Then by Theorem 11.31 we can arrange for v{xo, 9) to be any continuous function in 9. 
Therefore knowing Hy{xo) for all v solving (13.ip is equivalent to knowing 9 ■ Xu{xo, 9) 
for each 9 G Since this can be done for any xq G X, we can recover from 

knowledge of 9 ■ Vu{x, 9) for all {x, 9) E X x Then using the formula 

u{x,9) = u{x t+9,9) — f 9 ■ Xu{x + t9,9)dt (3-21) 

Jo 

and the fact that we know M|r+, we see that we can recover u{x, 9) at each point in 
X X Then using fll.ip . we can recover S. 

We now examine the stability of the reconstruction described above and thereby 
prove Theorem 11.11 


Proof of Theorem \1.1\ Suppose Si and S 2 are different sources. Let ui ^ and M 2 ,e be 
the solutions of fll.61) corresponding to Si and 5*2. Then since equation fl3.10l) applies 
to Ui^e and M 2 ,£, we have 


fax 


{ui^e — U 2 ,£)v n ■ 9 dxd9 


Jx 

+ 


(cr — CTg) {ui^s — U 2 ,£)v dx d9 + 


v{Si^f£ — S 2 ^e) dx d9 (3.22) 



'X 

If e = 0 we have 


'S"-i JX 

{ks — k){x, 9 ■ 9')v{x, 9){ui^s ~ U 2 ,£){x, 9') d9' d9 dx. 


/ / {uifl — U 2 fl)vn ■ 9 dx d9 = / / v{Si — S 2 ) dx d9. (3.23) 

Jax ’ ’ Jx 


Then fl3.22p becomes 


/ / {{ui^e - Uifl) - {U2,e - U2fl))v n ■ 9 dx d9 

Jax 

[ [ {a - as){ui^s-U2,£)vdxd9 + ( [ v{{Si^e - Si) - {S 2 ,£ - S 2 )) dx d9 


+ 



lx 


(/cg — k){x, 9 ■ 9')v{x, 9){ui^, 


U 2 ,e){x, d') d9' d9 dx. 


Sn-l 


(3.24) 
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Now note that 


'S'*-! JdX 


{ui^e — uifl)v n ■ 9 dx dO = / {ui^e — uifl)v n ■ 9 dx dO (3.25) 


= e 


'r+ 


n ■ 9 dx d9. 


(3.26) 


Therefore we can rewrite fl3.22p as 
(— Aljn n ■ 9 dxd9 


'S'*-! JdX 


cos(g • X + (p)a{ui — U2)v dx d9 + 



A Js**-! Js**-! 


/ V cos{q ■ X + ip){Si — S2) dx d9 
/S'*-! Ja Js**-! Jx 

cos{q ■ X + ^p)k{x, 9 ■ 9')v{x, 9){ui — U 2 ){x, 9') d9' d9 dx + 0{e). 

(3.27) 

Using the reasoning described above, we hnd that 

[ [ n ■ 9 dx d9 = — H ^^2 + 0{e). 

Js”-^ JdX 

It follows that for a hxed v solving (13.11) . we have 

II^Si “ ^S2lUh«"x{o,f},c'(aA))|k||Li(S'*-!,c(aA)) +0(s) > — 77 i;,2||li(R'*), (3.28) 

which means that 

\\Hv,i — Hy^ 2 \\L°°{x) < ||A|^ — As2lliMK'‘x{o,f},c(9A))lli'||L!(S'*-!,c(aA)) + 0{s). (3.29) 

Now by Theorem ll.3[ we can choose v snch that v{xo,9), as a function of 9, is an 
approximation of identity centered at 9o. Moreover, we have the estimate 

lklli,!(S'*-!,c(A)) < c \\ v { xo , 9 )\\ li ( s ^-^) < c, (3.30) 

where c only depends on X, a, and k. Then we obtain the estimate 

1^0 ■ V(mi — U 2 ){xo, 0o)| < c||^li “ ^S 2 llih®"x{o,f },c(9A)) + 0{e). (3.31) 

It follows by integration that 

r"^||-ui--U2||L°°(AxS'*-!) < c (^11A^^ - A5jlc(r+) + ||A|^ - A|J|ii(Rnx{o,|},c(SA)))+0(£) 

(3.32) 

and thus 

ells’! - 5*2IILOO(X) < IIA^^ - A5jlc(r+) + ||A|^ - A|J|Li(K'*x{o,f},c(9A)) + e(s), (3.33) 
as claimed. □ 

Note that if we are interested in the stability of reconstructing the source from the 
functional H^, we obtain the estimate 

eilS*! — S'2 ||l°°(a) < \\Hv,i — Sfj;,2||L°°(A) + HAg^ — A52l|c(r+)- (3.34) 
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In particular, the 0{e) error disappears, since this conies from the problem of recov¬ 
ering from A|. This is worth recording as a proposition, for comparison to the 
equivalent stability result in m- 

Proposition 3.2. If Si, S 2 are continuous functions on X x then for proper 

choice ofv, 

CII'S'i — S2\\l^(x) < \\Hv,i — IIv,2\\l°°{x) + — A52l|c(r+)- (3.35) 

where Hy^i, Hy ^2 are the functionals defined by Si and S 2 in terms of fl3.20p . 

4. Controllability 

The goal of this section is to prove Theorem ll.31 First we need the following lemma, 
which allows us to propagate solutions to larger and larger sets. 

Lemma 4.1. Suppose Xi d X 2 C X are concentric open balls and X 12 = X 2 \ Xi- 
Let 

T, = dX, X , (4.1) 

for i = 1, 2. Now suppose that ui is eontinuous on Xi, and solves 

9 ■ Xui + a{x)ui = I k{x,9,9')u{x,9') d9' (4.2) 

on Xi. Then for any p with 1 < p < 00 , there exists a solution u G L°°(X 2 )) 

to fl4.2p on X 2 such that u = Ui on Xi. Moreover 

\\u\\lp{S^-'^,L°°(X2)) < C\\Ui\\LP(S^-ifi{Xr)), (4.3) 

where C is a constant that depends only on a, k, and the diameters of Xi and X 2 . 

The proof of Lemma 14.11 requires a trace theorem. Similar trace theorems can be 
found in [H], [15], as well as [IB], [IS], and [1]. To state the trace theorem, we will 
once again use the notation introduced in the beginning of Section 2. 

Lemma 4.2. Ifu,TQU G L^{S"‘~^, L°°{X)), then for 1 < p < 00 

||n||LP(S"-i,L°°(9A)) < t\\Tqu\\lp(S’^-i,L°°(X)) + \\u\\lp{S"-'^,L°°{X))- (4.4) 

Proof. For almost every {x,6), we can write 

pT±{x,e) 

u{x ±T±{x, 9)9,9) = / Tou{x ±t9,9)dt + u{x,9). (4.5) 

Jo 

Now {x ± T±{x, 9)9, 9) E T± and dX x S^~^ = F+ U F_ up to a set of measure zero. 
Thus for almost every 9, 

|| n (-, 0 )|| L °°( aA ) < t\\Tou{-,9)\\loo(^x) + \\u{-,9)\\l<=o(x)- 

The lemma follows by taking norms in 9. 


Now we can prove Lemma [4.11 


(4.6) 

□ 











RADIATIVE TRANSPORT UMBLT 


11 


Proof of Lemma \4.1\ First, we will show that we can construct a function Ui 2 G 
L°°{X 12 )) which solves fl4.2p with the boundary condition Mi 2 |ri = Wilri- 
To do this, we will begin by picking a continuous function U 12 on X 12 x which 
satishes the boundary condition M? 2 lri = "Wilri and the estimate 

ll“l2lUp(S"-bC(Ai2)) < ||mi||lp(S’— i,C(9Ai))- (4.7) 

It follows from Lemma 14.21 that 

l|■*^12lk!’(S'*-^C(Al2)) < c\\Ui\\lp(^s^-i^c{Xi)), (4.8) 

for some constant c depending on a, k, and the possibly the diameter of Xi. 

Now we can dehne Ui 2 iteratively by solving the problem 


Tou^it^ 

= Au^2 oa ^12, 

(4.9) 

42+'|ri 

= "^^l Pi, 

(4.10) 

The solution to the general problem 



Tow 

= Q on X 12 , 

(4.11) 

^a|ri 

= hn. 

(4.12) 


of this form can be constructed explicitly by 

rr-(xfi) 

/ Q{x — t9,9)dt + h{x — t-{x,9)9,9), 


w{x, 9) 
w{x, 9) 
w{x, 9) 


'0 


rr+{x,e) 


Q{x + t9,9)dt + h{x + r+(x, 9)9,9), 


>0 


T-{X,0) 


Q{x — t9, 9)dt, otherwise. 


if x — T-{x, 9)9 G dXi 
if X + r+(x, 9)9 G dXi 

(4.13) 


Jo 

Here t± are dehned for (x, 9) G X 12 x with respect to the boundary of X 12 . 

Let A = max{r±(x, 9)\{x, 9) G X 12 }. When h = 0, we obtain from the expressions 
above that for any hxed 9, 


ll'^^(•l ^)IU°°(Vi2) < ^1 Q(-l ^)IU°°(Vi2)- 

(4.14) 

Therefore 


lh^l|LP(5"-i,L°°(Ai2)) A ^1 Q|LLP(S'‘-^L°°(Al2))• 

(4.15) 

NoW if Ui 2 are dehned iteratively as described above and are 

follows that 

defined by — 

^ 0 ^ 12 ^^ = ^v’l 2 on X 12 , 

(4.16) 

^^i"2+'|ri = 0. 

(4.17) 

Therefore 



(4.18) 


where 
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If A is sufficiently small, then a A < 1. Thus the converge geometrically to zero, 
and so converge in the Lp{S"’~^, L°°{Xi 2 )) norm to a function Ui 2 , where 

1 0 c 

||m12||lp(S"-i,L“(Ai 2)) < ^ ll'»12lUp(S"-hC°°(Ai2)) < ll^l ||lp(5’^-i,L°°(Ai))• 

(4.20) 

Since = Aui 2 , one can check that Tom ^2 converges in L°° as well, with a similar 

estimate applying. It follows that u\2 —t ni2 in IT = {n G L°°{Xi2)) : Tqu G 

L°°(Ai 2 ))}. Then ni 2 has trace ni 2 |ri = wilri by Lemma WfR so we can let 


f Mi(x, 9) if x G Xi 
\ Mi2(x, 9) if x G Ai 2 


(4.21) 


and check that u is a weak solution to fl4.2p . It follows that u must be the unique 
^p(^n-i, ^oo(j^^)) golution to fl4.2p on A 2 with ^^(5'”“^, L°°(r 2 ,_)) boundary condi¬ 
tions guaranteed by Theorem 12.11 Moreover, it is easy to see that u must satisfy the 
estimate fl4.3p since M 12 does. 

If A is not small enough to guarantee that a A < 1, we can repeat the above process. 
Note that if Xi and A 2 are concentric balls of diameter r and r + 6, respectively, then 

A = 2V2r(5 + (52. (4.22) 

Therefore if 1/a is sufficiently small, 

(5 = min I-,-I (4.23) 

\l0ar’l0aj ^ ^ 


would give a A < 1. Since the harmonic series diverges, it follows that only hnitely 
many repetitions of the above process are required to obtain a solution on a ball of 
any radius, satisfying the correct estimate. Note, however, that this implies that the 
constant C in the statement of Lemma 14.11 blows up as a or A go to inhnity. 

□ 


Proof of Theorem \1.3[ We begin by considering the case where ra < 1/2. Note that 
this occurs as long as X is sufficiently small. Now suppose Xq G X, and h{9) is 
continuous. We want to construct a solution v to the RTE with the property that 


v{xo,9) = h{9). To begin, consider the equation 

9-VvQ + a{x)vQ = f k{x,9,9')vo{x,9') d9', (4.24) 

Wo(2^,6')|r_ = h{9) (4.25) 

By Theorem 12.11 the above has a unique continuous solution vq with the estimate 

||'?^o||lp(S"-i,c(a)) < I — (4.26) 

Now let gi{9) = h{9) — Vq^Xq, 9). We can construct iteratively for j > 1 

9 ■ Xwj + a{x)wj = f k{x, 9,9')wj{x, 9') d9\ (4.27) 

Wj{x,9)\r_ = gj{9), (4.28) 
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where 

Now 


9j+i{0) = 9j{0) -Wj{xQ,e). 


9j{0) -Wj{xo,e) = 


r'T_(xo,0) 


'0 


{9 ■ Vwj){xo — t9,9) dt 


'•T_(XO,0) 


Then for a hxed 9, 
Therefore 


= / {Awj){xo — t9,9) dt. 

Jo 

\9j+iid)\ < T\\Awj{-,9)\\c(x)- 


(4.29) 

(4.30) 

(4.31) 

(4.32) 


||5'i+i|Up(S"-i) < 'r\\Awj{-,9)\\LP(s^-i,c(x)) < 'ra\\wj\\LP(^s’^-i^c{x))- (4.33) 

By Theorem 12.11 it follows that 

TQj 

\\9j+i\\Lp(s^-i) < — -!)• (4.34) 

X 1~ Qj 

Since we are assnming that ra < 1/2, we hnd that \\gj\\ip(^s’^-i) converges geomet¬ 
rically to 0. Using Theorem 12.11 again, we get||t(;j||LP( 5 n-i c'(x)) < ‘^\\9j\\Lp{s^-^), so 
follows that ||tr'j|lLP( 5 n-i c'(x)) converges geometrically to 0, and thns the snm 


v{x,9) = vo{x,9)+ ''^^Wj{x,9) (4.35) 

converges nniformly, by considering p = oo. The above yields a continnons fnnction 
V, which has the property that n(a;o, 9) = h{9), and which solves 


9-Vv + a{x)v= f k{x,9,9')v{x,9') d9' (4.36) 

Js”-^ 

with the bonndary condition 

v\t_ = h{9)+ J29j{G) (4.37) 

and the estimate 

lbllLP(S"-i,c(A)) < I — TOi + 1 _ 2ra (4.38) 

Now if TO > 1/2, then for xo G X, we can begin by taking a small ball Xq C X 

aronnd Xq. If the diameter of Xq is snfficiently small, then r^oO < 1/2, so by the 
above reasoning, there is a continnons fnnction v on Xq snch that v{xo,9) = h{9). 
Then by Lemma ITTl there is an L°°[X)) solntion to fl4.2p on some large ball 

containing X, which is eqnal to v on Xq. Restricting this to X hnishes the proof of 
Theorem 11.31 □ 


Remark 4.3. Note that since the constant C in the statement of Lemma \4.1\ blows up 
as a or X goes to infinity, it follows that the constant C in the statement of Theorem 
17.31 also blows up as a or the diameter of X become large. Thus the controllability of 
the RTE becomes worse as we approach the diffusion limit, which is to be expected. 
















14 


BAL, CHUNG, AND SCHOTLAND 


5. Two-Dimensional Constant Coefficient Case 


In this section we will consider a simplified version of the problem, where X C 
and the coefficients a and k are constant in x. In this case, the angular variable 6 now 
ranges over S^, and we can provide a different proof of Theorem 11.31 bv analyzing the 
RTE in terms of Fourier series in the angular variable. This simplihed approach may 
have applications to numerical experiments. It is possible that a similar approach 
with spherical harmonics works in three dimensions, but there the recursion relations 
are much more complicated, and it is not clear that the solutions converge. 

If we parametrize 6 & hy writing 6 = (cost, sint), for t G [0, 27r], then the RTE 
from Theorem 11.31 becomes 

p27t 

cos tdxV + sin tdyV + av = / k{t,t')v{x,t') dt'. (5.1) 

Jo 

We will prove the following version of Theorem 11.31 

Proposition 5.1. For each xq G and m G Z, there exists a solution v to fl5.1l) 
such that v{xq) = e*™* and ||t||l°o is bounded uniformly in m. 


Proof. We begin with the ansatz 

v{x,t) = ^T„(a;)e*”h (5.2) 

By the assumptions on the regularity of k, we can expand k{t,t') in terms of {e*”*} 
and {e*"* } as well. Crucially, the assumption that k obeys reciprocity implies that 
the expansion has the form 




(5.3) 


Therefore 


^271 

fc(t, t')v{x, t') dt' = knVn{x)e^"'^. 


Now we can rewrite fl5.ip as 


nEZ 


^•271 


e^^dzV + e '’^dzV P av = / k{t, t')v{x, t') dt', 


(5.4) 


(5.5) 


where dz = l/2{dx — idy) and dz = ll2{dx + dy). Substituting into the expansions 
for V and k and examining the e*""* term, gives the system of equations 


dz^n—l T dz^n+l i^kn a')Vn 


(5.6) 


Now suppose m G Z. Let Vn 


= 0 ioT n < m, Vm = 1, and 


Vn 


{n — m)\ 


(5.7) 


for n > m. The equations for e*"*, n < m are automatically satished, since each term 
in those equations is zero. For n> m, note that the dzVn-i is always zero. Therefore 
the equation for e*”* reads 


dzVn+l i^kn Cr'jVny 


( 5 . 8 ) 
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and one can check that it is satished by fl5.7j) . 

Lastly, the conditions on k imply that {/c„} are bounded, so there exists some 
constant C such that 

^n—m 

~ {n — m)\ 



(5.9) 


Therefore the sum 

v{x,t) = (5.10) 


converges and is bounded uniformly in m and v{x — Xq, t) a smooth solution to fl5.ip . 
such that v{xo,t) = □ 
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